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We evaluate the dynamic structure factor S{q,ui) of a one- dimensional quantum Hamiltonian 
with the inverse-square interaction (Calogero-Sutherland model). For a fixed small q, the structure 
factor differs from zero in a finite interval of frequencies of the width Suj oc q'^/m. At the borders of 
this interval S{q, to) exhibits power-law singularities with exponents depending on the interaction 
strength. The singularities are similar in origin to the well-known Fermi-edge singularity in the 
x-ray absorption spectra of metals. 
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Fermi liquid theory proved to be extremely successful 
in describing interacting fermions 1] . The low-energy ex- 
citations of a normal Fermi liquid (FL) are classified the 
same way as the excitations of a reference system - a non- 
interacting Fermi gas. A weak residual interaction leads 
to a finite decay rate 1/r of FL's quasiparticles. Formally, 
the rate can be defined as the width of the quasiparti- 
cle peak in the spectral function (imaginary part of a 
single-particle retarded Green function) A{q,u!). For FL, 
A(q, cj) at T = is a Lorentzian. Its width is easily eval- 
uated in the second order of perturbation theory, which 
yields 1/t cx uP'jf.p, where uj{q) is the quasiparticle en- 
ergy. 

It is well-known, however, that in one dimension (ID) 
even a weak interaction breaks down the FL description 
(see ^\ for recent reviews) . A guide to understanding the 
properties of interacting ID systems is provided by the 
Tomonaga-Luttinger model (TLM) |^, which plays the 
same role for the concept of the Luttinger liquid |7| as the 
Fermi gas does for FL. The TLM assumes a strictly lin- 
ear fermionic dispersion relation. With this assumption, 
the TLM Hamiltonian can be diagonalizcd exactly 0, 
no matter how strong the interactions are. The corre- 
sponding elementary excitations are bosons, quanta of 
the waves of fermionic density. These bosons do not in- 
teract, have an infinite lifetime, and propagate without 
dispersion, Lo{q) — uq (here u is the plasma velocity). 
Therefore, a measurable quantity, the dynamic structure 
factor (density-density correlation function) ^] 

S{q,uj)= fdxdte''-'^'~'i''\p{x,t)p{0,0)), (1) 



has an infinitely sharp peak, S'tlmI?,^) oc qS{uj — uq). 

Applications of TLM to the description of "real" ID 
fermions rely on the expansion of single-particle energies 
about the Fermi points ±pF 0,0, 



= ± vpk + k^/2m - 



PTPF, 



(2) 



where the upper/lower sign corresponds to the right/left 
movers (throughout this Letter we use units with h — 1). 
The linear in k term in Eq. (O is accounted for in TLM. 
The k^ term generates interaction between the bosons 



with the coupling constant oc 1/m which broadens 
the peak in S{q, uj). However, in contrast with FL theory, 
this broadening is inaccessible by perturbation theory. 

This can be seen by considering a special limit of LL, 
that of non- interacting fermions with spectrum In 
this case the structure factor differs from zero only if lo 
lies within a finite interval lo^ < lu < u^, where uj±{q) = 
vpq ^ q^/2m for q < 2pp. Within this interval S{q,uj) 
is constant, S = m/q. Thus, at q —> the structure 
factor indeed approaches the TLM form (with u = vp), 
but in a very peculiar fashion: the peak in S{q, lo) at 
a fixed q has a manifestly non-Lorentzian "rectangular" 
shape with the width Sto = q^jm. Even this simple result 
is non-perturbative in the bosonic representation: the 
first-order in 1/m contribution to the boson's self-energy 
vanishes, while the next one diverges on the mass shell 

Recently it was argued that in the presence of in- 
teractions between the fermions the shape of the peak 
in S'((7, lS) remains to be non-Lorentzian. Moreover, 
even a weak interaction transforms the discontinuities at 
LO = Li}± into power-law singularities. Here we approach 
the problem of finding S'(g, lo) from the perspective of the 
exactly solvable Calogero-Sutherland model (CSM). 

The CSM Hamiltonian reads [3 

N „2 

where Vix) is a periodic version of Xjx^ potential, 

A(A - l)/m 



V(x) 



{L/TT)'^sin\TTx/L) ■ 



(4) 



The excitations of CSM can be described in terms of 
quasiparticles and quasiholes Q . Quasiparticles are char- 
acterized by velocities v in the range \v\ > u, and an in- 
ertial mass m, the bare mass that enters Eq. Quasi- 
holes have velocities v in the range \v\ < u, and fractional 
inertial mass rh — m/\. The plasma velocity u is given 
by 8] 



TT\pf)/m, po = N/L. 



(5) 
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The momentum and energy (relative to the ground 
state) of an excited state of CSM characterized by a cer- 
tain set of velocities {vi, Vj} read Q 



P — mvi — mvj 

i J 



2 -2\ 



(6) 
(7) 



The result of the action of the local density operator 

L ^ N 

- ^ (5(x - Xi) - po 



pI= / dxe^'^'^pix), p{x 
Jo 



L 



on CSM's ground state |0) has a remarkable property (3: 
|0) is an eigenstate of CSM in which the velocities of 
all quasiparticles point in the same direction (positive 
for q > 0). This has a profound effect on the structure 
factor: S'(g, w) ^ only in a finite interval of frequencies, 
cj_(g) < a; < uj+{q), see Fig. QJa) (the upper bound uj+ 



PPe 

would be absent in a generic ID system [Tj). 



(a) 




(b) 




FIG. 1: (a) The structure factor S{q,u}) differs from zero in 
a finite interval of frequencies u!-<ij <ll!+. At the borders of 
this interval S{q,u!) exhibits power-law singularities, S (x — 
tj±|^*^"\ see Eqs. and The low-energy {lo- 0) 

sectors correspond to g = 2TTpoI, where I is an integer, (b) 
Dependence of S{q,iL!) on at a fixed q 7^ 2npoI and for a 
repulsive interaction (A > 1). 

From this point on, we consider the rational values of 
A only, A = r/s, where r and s are co-primes. In this 
case, the state \R) — p^^Q |0) has exactly s right -moving 



quasiparticles and exactly r quasiholes 0,0) This is 
the simplest possible excitation that conserves the total 
inertial mass, 



sm — rm — for A r / s. 



(8) 



The bounds LiJ± are easily found from the momentum and 
energy conservation, 



Pi 



LU^ El 



R)- 



(9) 



Indeed, it follows from Eqs. that 

s r 

q = m{vi ~ u) + fh{u — Vj), (10) 



J = l 



2{llJ — uq) = m{vi 



1=1 



J2rhi^~v,f. (11) 



Since Vi — u>0 and u — vj > 0, Eq. H10() implies that for 
a given q < 27rpo the velocities vary in the range 



u < Vi < Vq, vq — u + q/m, 
Vq < Vj < u, Vq — u — q/fh. 



(12) 



The upper bound w = is reached when the velocity 
of one of the quasiparticles approaches Vq while all the 
remaining quasiparticles/holes have velocities close to the 
plasma velocity u. Eq. Hll() then gives 



uq^^{vQ-u) . 



(13) 



Similarly, the lower bound corresponds to the inter- 
mediate state \B) in which one of the quasiholes has ve- 
locity close to Wo, while the velocities of all the remaining 
quasiparticles/holes are close to u, so that 



a;_ — uq ■ 



(u-vof 



q 

2m 



-A- 



2 ' ' 2m 2m 
The width of the region where 5 7^ is then given by ^ 



(14) 



5uj 



(A + 1) 



2m 



(15) 



Eq. H14|l is valid as long as w_ > 0, i.e. for q < 2-KpQ. 
At larger q the function Ld^{q) is given by Eq. H14|) with 
the replacement q q — iirp^I , where / is the integer 
part of qjl-KpQ. The corresponding intermediate state 
|_R) has / quasiholes with velocities approaching — u, one 
quasihole with velocity near vq (given by Eq. H12I) with 
the replacement q ^ q — 27tpqI), and all the remaining 
quasiparticles/holes moving with the plasma velocity u. 

We turn now to the evaluation of the structure fac- 
tor. In the thermodynamic limit (TV 00, po = const) 
Eq. Q can be rewritten as 



S{q,uj) ^ W Fs,rS{q - P^r)) S(uj - E\ 



(16) 



The form-factor Fs^r here is given by 



Fg r OC 



Hi 



\Vi - Vi 



|2A 



, l-A 



(17) 

(note that Y\dvidvjFs^r is dimensionless) . This expres- 
sion was conjectured in |0| based on the results of 
Ref. for A = 1/2, 2 the conjecture was proved in 
using properties of Jack polynomials. In writing H17|l . 
we omitted A-dependent numerical coefficient Qj. 

For simplicity, we concentrate on the most interesting 
limit of small q <C Trpo- In this limit vq — u,u — vq u, 



see Eq. (|12ll , and one can approximate 



2u{vi 



m), u 



2u{u — Vj) in Eq. H17|l . In view of the 
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restriction H12() on the velocities, it is convenient to switch 
in Eq. (|16|l to the new integration variables 



U — V4 



(18) 



Vq ~ U U — Vq 

which vary between and 1. In terms of these variables 

S{q,uj) 



So{q) 



^1 

J + (19) 

LU — uq 



Ea;?-AEst-(A + l 



Slo 



where So{q) — m/q is the structure factor for noninter- 
acting fermions (A = 1), and the form-factor is 



n 



|2A 



i<ii' 



n 



3<j' 



(20) 



3 3 



Evaluation of Eq. H19|l simplifies considerably when 
uj u!±. Consider, for example, the limit w — u- Slo. 
In this limit the non-vanishing contributions to the inte- 
gral in (Ull come from r sectors of the (r-l-s)-dimensional 
integration space in which one of Xj, say Xfe, is close to 1, 
while all the other variables approach 0. To lowest order 
in Xi, Xj^k, 1 — 5fe *C 1, the second (5-function in Eq. I|19(l 
can be replaced by 



Slo 



(5(...)« S[2\il-Xk)-iX + l) 
and the form-factor F^^r becomes 

Fs,r(^{x^},{xj}'j W X-^'Fs,r-l(^{x^},{xJ}J 



{Fs.r-1 is given by (|20|l with j = k excluded). With 
these approximations, the integrations in H19I) are easily 
carried out and yield 



Soiq) 



oc 



Slj 



l/A-l 



< LO - LO- 't: Slo. (21) 



Similarly, at 0;+ — lo <C Slo, the non- vanishing contri- 
butions to the integral in Eq. (|19|l come from s sectors 
where 1 — Xk,Xi^k,Xj ^ 1, and one finds 



S{q,uj) 
Soiq) 



LO^— LO 



Slo 



A-l 



< L0+ ~ LO <t: Slo. (22) 



It should be emphasized that Eqs. (|21|l and H22|l with 
q-independent exponents are valid for all q ^ 2'kpqI, in- 
cluding q > npo. In addition to Lo±{q), a smooth depen- 
dence on q enters Eqs. (1^ and via the (omitted) 
prefactors; these w-independent prefactors tend to 1 in 
the limit A ^ 1. According to Eqs. (|2lJ and (|221), the 
structure factor diverges at uj co- {uj+) for a repulsive 
(attractive) interaction, see Fig. ^b); the divergencies 
are integrable for all physically meaningful (i.e. positive) 



values of A Note that, formally, Eq. (|22(l can be ob- 
tained by replacing \lo — w_| ^ \lo — lo^\ and A — > 1/A in 
Eq. (|21|l . This correspondence is a manifestation of the 
particle ^ hole, A^l/A duality of CSM (li |. 

The power-law singularities in Eqs. H21II and (|22|l are 
similar to the familiar edge singularities in the x-ray ab- 
sorption spectra in metals |l7|. Consider, for example, 
the limit lo — » lo+. In this case the action of the operator 
p}j on the ground state creates a quasiparticle at k fa q 
with energy « loj^. This high-energy quasiparticle in- 
teracts with the Fermi sea, resulting in the excitation of 
multiple low-energy particle-hole pairs, just like the core 
hole does in the conventional x-ray edge singularity. The 
proliferation of the low-energy particle- hole pairs leads to 
the singularity in the response function 0, E| . 

This analogy suggests that, just like in the case of the 
conventional edge singularity, the functional form of the 
dependence of the structure factor on lo — can be cap- 
tured by replacing the original model with the properly 
chosen effective Hamiltonian. For lo — *■ a;+ and q irpo 
it is sufficient to include in the effective Hamiltonian only 
the right-moving single-particle states within two very 
narrow stripes (subbands) of momenta near k = q and 
fc = to allow for both the creation of the high-energy 
quasiparticle and the low-energy particle-hole pairs. (Re- 
call that for q < irpo the velocities of all quasiparti- 
cles/holes in the state |0) are positive, see Eq. (I12|l : 
this would not be the case for a generic interaction 7]). 

Upon introducing 



|fe|<fco 



fc-g|<fco 



■ ll'k, 



where ipk annihilates a right-moving particle with mo- 
mentum fc (fc = at the Fermi level) and fco ^ 9 is a 
cutoff, the effective Hamiltonian can be written in the 
coordinate representation. 



H+ = dxTp. 



dx ip] 



U I dx Pd{x)pr{x). 



(23) 



Here pr^d = ■ d^r d ■' where the colons denote the nor- 
mal ordering. In Eq. (|23|l the nonlinearity of spec- 
trum (O is encoded in the mismatch of the velocities 
vq ~ u ~ q/m, see Eq. ()12|l. The inter-subband interac- 
tion constant U is set by the requirement that the two- 
particle scattering phase shift Q = U / (vq — u) for Eq. H23fl 
is equal to that for CSM, 6 = (A - 1)tt 8]. This gives 

U ={vo- u)e = (A - l)nq/m, (24) 

which for |A - 1| < 1 coincides with Vq -Vq fHK^. 
In terms of Eq. H23|l , the structure factor is given by 

S{q,Lo)^ fdxdte''^'{bix,t)b\0,0)), 6^ = VdV'r- (25) 
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Note that the total number of d-particles A''^ = Jdx Pd{x) 
commutes with and that the entire d-subband hes 
above the Fermi level. Hence, as far as the evaluation of 
Eq. is concerned, can be further simplified by 
replacing ipdix) — *■ V'ipdix)V, where P is a projector onto 
states with Nd = 0, 1. Obviously, the projected operators 
satisfy pd{x)'ipd[y) = 0, and i>dix)pd{y) = S{x - y)ipd[x), 
which implies that \_pd{x)^ pd{y)\ = 0. 

We now bosonize the f/'r-field according to 

i!r{x) = ^/h) e''P'^''\ [p{x),ip{y)\ = i7rsgn(x- y), 

and apply a unitary transformation 7] with generator 
W = {Q /2TT)j dx pd{x) dxip- The transformed Hamilto- 
nian reads 

H+ = e'^H+e-'^ = i^o + SH, (26) 

where 

Ho = ^J dxid^Lpf + j dx^lj\{Lu+- ivodx)^:^ (27) 

and SH — {vq — u){Q'^/ATT)Jdxp^{x). It is easy to see 
that a state with a single d-particle is an eigenstate of 
SH, dH^jl{x)\0) cx {ko/q)6uj{X-l)^^jl{x)\0). Thus, when 
acting in the subspace with Nd = 0, 1, the second term in 
Eq. (Uni results merely in a correction to uj^ in Hq , which 
for ko/q 1 can be safely neglected, i.e. « Hq. 

The same unitary transformation applied to the oper- 
ator fet in Eq. if^ yields 

V{x) = e'^^b^x) e-'^ = ^^o^^lix) e<i+e/2-)^W. 

Evaluation of the correlation function (|25|l with 
quadratic Hamiltonian Hq is now straightforward. The 
structure factor vanishes identically at w > while at 
Lu < uj-f- it is given by Eq. H22|l . Thus, the outlined simpli- 
fied description indeed reproduces the exact result 
with logarithmic accuracy. (The cutoff fco would enter 
Eq. 122|l via a factor g/fco in the square brackets.) Note 
that the exponent in Eq. (|22|l is independent of q. This 
independence is a direct consequence of the fact that the 
phase shift & — const for inverse-square interaction. 

Similar reasoning can be applied to the calculation of 
S{q,uj) at — > uj-. In this case the c?-subband lies near 
k — —q well below the Fermi level and carries at most 
a single hole 0- After the particle- hole transformation 

d ^ ^1 d corresponding effective Hamiltonian i7_ 
takes the form of Eq. I|23|l with replacements uj^ w_ 
and vq —> vq. Evaluation of S{q,uj) [which is again given 
by Eq. if^ ] proceeds similar to above and yields Eq. (|^ . 

To conclude, in this Letter we evaluated the dynamic 
structure factor S{q,uj) of the Calogero-Sutherland 
model. Besides being of a fundamental interest for it's 
own sake, the detailed knowledge of the structure fac- 
tor for interacting fermions with nonlinear dispersion is 



important for the description of a variety of effects asso- 
ciated with the particle-hole asymmetry. 

We found that S{q, lo) differs from zero in a finite inter- 
val of frequencies. At the borders of this interval S{q, lu) 
exhibits power-law singularities, analogous to the edge 
singularities in the x-ray absorption spectra of metals. 
Exploiting this analogy, we showed that the exact re- 
sults (jnj and (|23l can be reproduced with logarithmic 
accuracy by replacing the original model with simple ef- 
fective Hamiltonians. Remarkably, the analogy with the 
x-ray singularity, previously established for weak interac- 
tion only 0], is useful even when interactions are strong. 
Moreover, similar ideas can be applied to the evaluation 
of single-particle correlation functions ^2^ . 
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